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Q ■ Abstract: We study the impact of flavour in "soft leptogenesis" (leptogenesis induced 
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by soft supersymmetry breaking terms). We address the question of how flavour effects 
can affect the region of parameters in which successful soft leptogenesis induced by CP 
violation in the right-handed sneutrino mixing is possible. We find that for decays which 
occur in the intermediate to strong washout regimes for all flavours, the produced total 
5^ I B — L asymmetry can be up to a factor 0(30) larger than the one predicted with flavour 

effects being neglected. This enhancement, permits slightly larger values of the required 
lepton violating soft bilinear term. 
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1. Introduction 

The discovery of neutrino oscihations makes leptogenesis a very attractive solution to the 
baryon asymmetry problem [1,2]. In the standard framework it is usually assumed that the 
tiny neutrino masses are generated via the (type I) seesaw mechanism [3] and thus the new 
singlet neutral leptons with heavy (lepton number violating) Majorana masses can produce 
dynamically a lepton asymmetry through out of equilibrium decay. Eventually, this lepton 
asymmetry is partially converted into a baryon asymmetry due to fast sphaleron processes. 

For a hierarchical spectrum of right-handed neutrinos, successful leptogenesis requires 
generically quite heavy singlet neutrino masses [4], of order M > 2.4(0.4) x 10^ GeV for 
vanishing (thermal) initial neutrino densities [4, 5] , although flavour effects [6-9] and/or 
extended scenarios [10, 11] may affect this limit. The stability of the hierarchy between 
this new scale and the electroweak one is natural in low-energy supersymmetry, but in 
the supersymmetric seesaw scenario there is some conflict between the gravitino bound on 
the reheat temperature and the thermal production of right-handed neutrinos [12]. This 
is so because in a high temperature plasma, gravitinos are copiously produced, and their 
late decay could modify the light nuclei abundances, contrary to observation. This sets an 
upper bound on the reheat temperature after inflation, Trh < 10^^^*^ GeV, which may be 
too low for the right-handed neutrinos to be thermally produced. 

Once supersymmetry has been introduced, leptogenesis is induced also in singlet sneu- 
trino decays. If supersymmetry is not broken, the order of magnitude of the asymmetry and 
the basic mechanism are the same as in the non-supersymmetric case. However, as shown 
in Refs. [13-15], supersymmetry-breaking terms can play an important role in the lepton 
asymmetry generated in sneutrino decays because they induce effects which are essentially 
different from the neutrino ones. In brief, soft supersymmetry-breaking terms involving 
the singlet sneutrinos remove the mass degeneracy between the two real sneutrino states of 
a single neutrino generation, and provide new sources of lepton number and CP violation. 



As a consequence, the mixing between the two sneutrino states generates a CP asymme- 
try in the decay, which can be sizable for a certain range of parameters. In particular 
the asymmetry is large for a right-handed neutrino mass scale relatively low, in the range 
10^ — 10* GeV, well below the reheat temperature limits, what solves the cosmological 
gravitino problem. Moreover, contrary to the traditional leptogenesis scenario, where at 
least two generations of right-handed neutrinos are required to generate a CP asymmetry 
in neutrino/sneutrino decays, in this new mechanism for leptogenesis the CP asymmetry 
in sneutrino decays is present even if a single generation is considered. This scenario has 
been termed "soft leptogenesis", since the soft terms and not flavour physics provide the 
necessary mass splitting and CP-violating phase. 

In general, soft leptogenesis induced by CP violation in mixing as discussed above 
has the drawback that in order to generate enough asymmetry the lepton-violating soft 
bilinear coupling has to be unconventionally small [13, 14]. Considering the possibility of 
CP violation also in decay and in the interference of mixing and decay of the sneutrinos [15], 
as well as extended scenarios [16, 17], may alleviate this problem. 

In Refs. [13-15] soft leptogenesis was addressed within the 'one-flavour' approximation. 
This one-flavour approximation is rigorously correct only when the interactions mediated 
by charged lepton Yukawa couplings are out of equilibrium. This is not the case in soft 
leptogenesis since, as mentioned above, successful leptogenesis in this scenario requires a 
relatively low right-handed neutrino mass scale. Thus the characteristic T is such that the 
rates of processes mediated by the r and /i Yukawa couplings are not negligible implying 
that the effects of lepton flavours have to be taken into account. 

The impact of flavour in thermal leptogenesis in the context of the standard see- 
saw leptogenesis has been recently investigated in much detail. [6, 7, 9, 11, 18, 19, 21, 23, 
24]. The relevant Boltzmann Equations (BE) including flavour effects associated to the 
charged lepton Yukawa couplings were first introduced in Ref. [18]. Additional flavour 
effects associated to the light-to-heavy neutrino Yukawa couplings which are particularly 
relevant for the case of see-saw resonant leptogenesis were discussed in Ref. [19]. In Ref. [6, 
7,11] it was further analyzed how flavour effects can significantly affect the result for the 
final baryon asymmetry. 

In this work we study the impact of flavour in soft leptogenesis. We address the 
question of how flavour effects can affect the region of parameters in which successful 
leptogenesis induced by CP violation in the right-handed sneutrino mixing is possible, 
and in particular their impact on the required value of the lepton-violating soft bilinear 
coupling. The outline of the paper is as follows. Section 2 revisits the soft leptogenesis 
scenario with CP violation in the mixing and we present the relevant BE describing the 
production of the lepton asymmetry in this scenario without including flavour effects. In 
Sec. 3 we discuss the way to include flavour-dependent processes associated with the lepton 
Yukawa couplings in this scenario. Finally in Sec. 4 we present our quantitative results. 



2. Unflavoured Soft Leptogenesis 

The super symmetric see-saw model could be described by the superpotential: 

W = ^M.jNiNj + Yi^eapN^L'^H^, (2.1) 

where i,j = 1,2,3 are flavour indices and Ni, Li, H are the chiral superfields for the RH 
neutrinos, the left-handed (LH) lepton doublets and the Higgs doublets with e^p = —e^a 
and ei2 = +1- The corresponding soft breaking terms involving the RH sneutrinos Ni are 
given by: 

Csoft = -m%N*Nj - ( AijYije^pNiqh^ + ^Bi^M^jN^N, + h.c.) , (2.2) 



where if = i VijiJ 1 and h = [h^ , h^j are the slepton and up-type Higgs doublets. 

The Lagrangian for interaction terms involving RH sneutrinos Ni and RH neutrinos 
Ni in 4-component spinors is given by: 

Qnt = -Yije^t3 [MiNfe^h" + PPd'-Ni + PPiN^t + ANith'') + h.c. (2.3) 



where £j = (z^j, £~) ,h^= (h^,h^] are the lepton and fermionic partner of h and Pl,r are 
the left or right projection operator. 

The sneutrino and antisneutrino states mix with mass eigenvectors 

iV+, = i=(e**/27v, + e-**/2iV;), 

iV_, = ^(e^*/2iV, - e-'^^/^N*), (2.4) 

V2 

where ^ = arg(BM) and with mass eigenvalues 

Mii^ = Mf, + ml±\BuMii\. (2.5) 

From (2.3) and (2.4), we can write down the Lagrangian in the mass basis as 



+iN-i 



h^PLi'^ + (Aj - Mi)qh'^'\ } + PPiNi^ + h.c. (2.6) 



In what follows, we will consider a single generation of A^ and N which we label as 1. 
We also assume proportionality of soft trilinear terms and drop the flavour indices for the 
coefficients A and B. As discussed in Refs. [13,14], in this case, after superfield rotations 
the Lagrangians (2.1) and (2.2) have a unique independent physical CP violating phase: 

<j) = s.Tg{AB*) (2.7) 

which we chose to assign to A. 

Neglecting supersymmetry breaking effects in the right sneutrino masses and in the 
vertex, the total singlet sneutrino decay width is given by 



Y,\M\\Y, 



,2 
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^N, = r^_ - r^ = ^-^- . (2.8) 



2.1 The CP asymmetry 

As discussed in Ref. [14], when T ^ AM± = Mj^ — M_, the two singlet sneutrino states 
are not well-separated particles. In this case, the result for the asymmetry depends on 
how the initial state is prepared. In what follows we will assume that the sneutrinos are in 
a thermal bath with a thermalization time T^^ shorter than the typical oscillation times, 
AAf_^r , therefore coherence is lost and it is appropriate to compute the CP asymmetry in 
terms of the mass eigenstates Eq.(2.4). 

As we will see below, the CP asymmetry produced in the decay of the state Ni=± 
which enters into the BE is given by: 

U = ^^-= , (2.9) 

where a^ = s^, /^ with s^ = ikh and /^ = i^h and we denote by 7 the thermal averaged 
rates. For convenience we also define the fermionic and scalar CP asymmetries in the decay 
of each Ni as 



(2.10) 



^/- = vA,.-. .c. ..„ ,.^. ..-. F7^ • (2-11) 



Notice that e^. and e/- are defined in terms of decay amplitudes, without the phase-space 
factors which, as we will see, are crucial to obtain a non- vanishing CP asymmetry [13,14]. 
Neglecting supersymmetry breaking in vertices, the total asymmetry e^ generated in 
the decay of the singlet sneutrino Ni can then be written as: 

e, = '^'^^' + '/'^/^ , (2.12) 

where Cs^jCf- are the phase-space factors of the scalar and fermionic channels, respectively. 
We compute the CP asymmetry following the effective field theory approach described 
in [28], which takes into account the CP violation due to mixing of nearly degenerate 
states by using resumed propagators for unstable (mass eigenstate) particles. The decay 
amplitude Aif of the unstable external state Ni defined in Eq. (2.4) into a final state a is 
described by a superposition of amplitudes with stable final states: 

j-^i I J JJ 
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where A^" are the tree level decay amplitudes and Hij are the absorptive parts of the 
two-point functions for i,j = it. The amplitude for the decay into the conjugate final state 
is obtained from (2.13) by the replacement A4" -^ A4"*. 

Neglecting supersymmetry breaking in vertices and keeping only the lowest order con- 
tribution in the soft terms we find the known result [13, 14] 



-e/+ 



-e/- 



Im^ AFB 



(2.14) 



As long as we neglect the zero temperature lepton and slepton masses and small Yukawa 
couplings, the phase-space factors of the final states are flavour independent and they are 
the same for i = it. After including finite temperature effects in the approximation of 
decay at rest of the Nj- they are given by: 



Cf^T) = CfjT) ^ Cf{T) = (1 - X, - x^)A(l,x,,x^) [l - f^'^] |l - f. 
Cs^{T) = Cs_{T) = Cs{T) = X{l,Xh,x^) [1 + /f ] 



fcq 
'h 



l + /~ 



where 



feq _ ^ 

'm exp[ii;. ./r]-i 






eq 

hi 



1 



exp[ii;^^,/r] + 1 

are the Boltzmann-Einstein and Fermi-Dirac equilibrium distributions, respectively, 

^(.~h~ ^\^^ ^Lh~ ^hl)^ ^1 



'hi 



2(1 + Xf^^l Xl^jJ 



niajTY 



A(l, X, y) = a/(1 + X - y)2 - 4x, Xa = — ^g^- 
The thermal masses for the relevant supersymmetric degrees of freedom are [29]: 
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2.15) 
2.16) 

2.17) 
2.18) 

and 

2.19) 
2.20) 



ml{T) = 2m|(T) 
m){T) = 2mj{T) 



l,l + \,l + l>?AT^ 



1 



92 + ^9y T' ■ 



2.21) 
2.22) 



Here §2 and gy are gauge couplings and At is the top Yukawa, renormalized at the appro- 
priate high-energy scale. 

As we will see in Sec. 2.2 the contribution to the relevant BE for the lepton number 
scalar and fermion asymmetries can be factorized respectively as: 



^7(iV±^Sfc)-7(A^±^Sfc) 



UT) 



Cs{T) 



J^7(iV±^afc)+7(iV±^afc) Cs{T)+Cf{T) 



(2.23) 



ah.,k 



and 

rrp\ — _^ _ - ^f\^) (2 24) 

The total CP asymmetry generated in the decay of any of the sneutrino N± is then: 

^(^) = ' 'irir'm " ' ^^^(^) • (2.25) 

Cs{T)+Cf{T) 

In this derivation we have neglected thermal corrections to the CP asymmetry from the 
loops, i.e., we have computed the imaginary part of the one- loop graphs using Cutkosky 
cutting rules at T = 0. These corrections are the same for scalar and fermionic decay 
channels, since only bosonic loops contribute to the wave- function diagrams in both cases, 
so they are not expected to introduce significant changes. 

2.2 The Boltzmann Equations 

We next write the relevant classical BE describing the decay, inverse decay and scattering 
processes involving the sneutrino states. 

As mentioned above we assume that the sneutrinos are in a thermal bath with a ther- 
malization time shorter than the oscillation time. Under this assumption the initial states 
can be taken as being the mass eigenstates in Eq. (2.4) and we write the corresponding 
equations for those states and the scalar and fermion lepton numbers. The CP fermionic 
and scalar asymmetries for each Ni defined at T = are those given in Eq. (2.14). 

The BE describing the evolution of the number density of particles in the plasma are: 



dnx 
dt 



+ 3Hnx = Y.Af^-[fiU...{l±fx){l±fj)...W{lm---^Xj..: 

j,l,m 

- fxfj . . . (1 ± //)(1 ± /™) . . . W{Xj . . . ^ /m . . . )] 



where, 






d?Px f d^Pj f d^pi f d^p. 



{2ttY2Ex J {27rf2Ej "J {27r)^2Ei J {2tt)^2E, 



and W{lm ■ ■ ■ ^ Xj . . . ) is the squared transition amplitude summed over initial and final 
spins. In what follows we will use the notation of Ref. [25]. We we will assume that the 
Higgs and higgsino fields are in thermal equilibrium with distributions given in Eqs. (2.17) 
and (2.18) respectively. Strictly speaking this implies that we are not including all the 
effects associated with spectator processes [26,27]. For the leptons and sleptons we assume 
that they are in kinetic equilibrium and we account for their asymmetries by introducing 
a chemical potential for the leptons, fi£, and sleptons, fif 

f^ = e{s,-w)/T + i' f^ = g(£,-M,-)/T _ 1 ' (2.26) 



and the corresponding ones for the antiparticles with the exchange ni —> —fj,£ and fi^ — > —fi^ 
respectively. Furthermore in order to ehminate the dependence in the expansion of the 
Universe we write the equations in terms of the abundances Yx, where Yx = nx/s. Also 
for convenience we use the variable z = M/T. 

We are interested in the evolution of sneutrinos Y^ , and the fermionic Yj^, and scalar 

Moreover, in order to 



yi^: 



Yl^ lepton numbers, defined as Yl^ = {Ye — Y^), Yl^ = (Y^ 

account for all the AL = 1 terms we also need to consider the evolution of the right-handed 

neutrino Yat. 

Neglecting supersymmetry breaking effects in the right sneutrino masses and in the ver- 
tices, all the amplitudes for Nj^ and N^ are equal as well as their corresponding equilibrium 

fl^ = f^S. So we can define a unique BE for Yr-r =Y^j +l"rr . 
■'TV- ■'Af ^ JV^ot + 



eq 



number densities, /- 

Thus, in total, in this unflavour case, we have a set of four BE. 

The derivation of the factorization of the relevant CP asymmetries including the ther- 
mal effects is somehow lengthy but straight forward. In particular one has to use that at 
0{e) we can neglect the difference between f^ and f^~ in the definitions of the ther- 
mal average widths (see for example Ref. [16]). Many of the terms in the equations are 
equivalent to the ones given for example in Ref. [31]. * 

Altogether we find : 



sHz 



dYN 
dz 



sHz- 



d^N. 



tot 



dz 



sHz — —^ 

dz 



(^-1 



/Yc-r 



IN + 47f ) + 47« + 47p) + 27f ) + 47^^ 



(2.27) 



N, 



y£9 

N 



^ - 2 (7^ + jf + 37.. + 27f ) + 27f ) + 27r) + 7^ ^ + 27^^ 



-7jv 

7jv 

_Yl 
Y 



YL,ef{T)+YL,es{T) 



^}{T) 






(2.28) 



tot 



N 



eq 



-;1N + 



^tot .,(5) 



y''iN 



.(6) 



„(7) 



Y, 



^ .,(3) 



y£9 

N 



7r + 27r+27r^+T7^7r+27i^ 



.(4) 



^N 



, Yl, - Yl^ 

+ peg 7MSSM' 



(2.29) 



sHz 



dYL, 
dz 




'^N. 



tot 



ye. 



Y-' 

N 



,.(3) 



-^IN +7^+2 yeg 

N 



V y''^N 

l^^tot (8) 



„(9) 



Yn (0) 



„(i) 



7r + 27r+2-^7r + 27r+27t^ 



„(2) 



^N 







7MSSM 



(2.30) 



'However, some care has to be taken as the Eqs. in Ref. [31] are given in the weak basis for the N while 
we give here the corresponding equations in the mass basis. 



In the equations above Yc'^ = j^ and K?''(r > M) = 90C(3)/(47r^5*), where g* is the 
total number of entropic degrees of freedom, (7* = 228.75 in the MSSM. 
The different 7's are the thermal widths for the following processes: 

722 = -f{N±i^uq) = -f{N±q* ^ tu) = j{N±u* ^ tq), 

7p) = 7(iV±£ ^ qu) = -/{N±i ^ qu) , 

7^ ) = ^(iV±u ^ l"g) = 7(iV±g* ^ lu) , 

7^(^) = 7(iV±^ ^ l"n) = 7(iV±n ^ I"*?), 

^f ) = ^{Nj* ^qu), 

7^ ) = 7(iV±g ^ ^) = 7(iV±n ^ ^^), 

7^ = 7(iV ^ ^/i) + 7(7V ^ th), 

7^ ) = ^(ivl ^ ^u) = ^(iV£ ^ qu), 

7W = 7(iV^ ^ tu) = 7(iV ^ rg) , 

7p) = 7(iVu* ^ tq) = -f{Nq* ^ tu) , 

7^ ) =^(^Ni^ qu) , 

7^^) = 7(iV ^ £g) = -f{Nq ^ lu) , (2.31) 

where in all cases a sum over the CP conjugate final states is implicit. 

We have included in Eqs.( 2.27-2.30) the N± and N decay and inverse decay processes 
as well as all the AL = 1 scattering processes induced by the top Yukawa coupling. We 
ignore AL = 1 scattering involving gauge bosons. We have accounted for the dominant 
CP asymmetry in the mixing as generated by the thermal effects in the N± two body 
decays but we have not included the possible CP violating effects induced by mixing in its 
three body decays or in its scattering processes. AL = 2 processes involving the on-shell 
exchange of N or N± are already accounted for by the decay and inverse decay processes. 
The AL = 2 off-shell scattering processes involving the pole-subtracted s-channel and the 
u and t-channel, as well as the the L conserving processes from N and N pair creation and 
annihilation have not been included. The reaction rates for these processes are quartic in 
the Yukawa couplings, ie they involve factors {YY'^)'^ , and therefore can be safely neglected 
as long as the Yukawa couplings are much smaller than one, as it is the case. 

The explicit expressions for the 7's in Eq. (2.31) can be found, for example, in [31] 
for the case of Boltzmann-Maxwell distribution functions and neglecting Pauli-blo eking 
and stimulated emission as well as the relative motion of the particles with respect to the 
plasma \ With these approximations, for example: 

eq-r. K.l{z) eq-r. K-l{z) /o QO\ 



^Neglecting supersymmetry breaking effects in the right sneutrino masses and in the vertices, it can be 
shown that the thermal widths for the sneutrino mass eigenstates and weak eigenstates are the same 



where /Ci^2(-2) are the modified Bessel function of the second kind of order 1 and 2 and 
r^v = r^ are the zero temperature widths Eq. (2.8). In our calculation we keep the thermal 
masses and statistical factors on the CP asymmetries but we neglect them in the rest of 
the thermal widths, with the exception of the Higgs mass the in the AL = 1 processes 
involving a Higgs boson exchange in the t-channel. 

iMSSM represent processes which transform leptons into scalar leptons and vice versa 
(for example [e + e •s-^ e + e]). The rates for these reactions are larger than the ones 
in Eq. (2.31) because they do not involve the Yukawa couplings Yij. Consequently they 
enforce that Y^j f» Y^^ . 

For Yl = Yl^ we can combine the BE for Y^, and Y^^ by defining 

Yl,^^^Y,^.+Y,^, (2.33) 



which obeys the BE: 
dYL, 



sHz—^ 
dz 



e{T) 



^tot _9 1 _ -^^tot 
Y^^ ] ^Yc" 



Tat 



"^ <'lN + ^7f ^ + 2,f + 2,?^ + ^,f ) + 2,f^^ 



2Y^^ I 4 Y~^ ' ^ It I ^ It < Y' 



' N 



eq 

N 




"^ ''iW + ^f + f^^?' + 27P' + 2|^/> + 2,1" + 2,f' ) 

N ^ / 

722. (2.34) 

Also the second line in Eq. (2.28) can be written as —^^e{T) „y° . So in total we are left 

with three BE for Ym , 1"^; , and Yr , , . 

^tot tot 

The final amount of B — L asymmetry generated by the decay of the singlet sneutrino 

states assuming no pre-existing asymmetry can be parameterized as: 

Yb-Uz ^ oo) = -Yl,Az ^ oo) = -2rjeY'^\T » M) (2.35) 

where e is given in Eq.(2.14) -t. 

r/ is a dilution factor which takes into account the possible inefficiency in the production 
of the singlet sneutrinos, the erasure of the generated asymmetry by L-violating scattering 
processes and the temperature dependence of the CP asymmetry and it is obtained by 
solving the array of BE above. Within our approximations for the thermal widths, r/ 
depends on the values of the Yukawa couplings (yy^)ii and the heavy mass M, with the 
dominant dependence arising in the combination 

{YY^)iivl = m,ffM (2.36) 



'"The factor 2 in Eq. (2.35) arises from the fact that there are two right-handed sneutrino states while 
we have defined V^ for one degree of freedom. Defined this way, r) has the standard normalizaiton r; ^ 1 
for perfect out of equilibrium decay. 
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Figure I: Efficiency factor |ry| as a function of rrieff for M = 10^ GeV and tan/3 = 30. Tlie 
two curves correspond to vanishing initial N abundance (solid black curve) and thermal initial N 
abundance, (dashed red curve). 



where Vu is the vacuum expectation value of the up- type Higgs doublet, Vu = v sin/3 (f=174 
GeV) . There is a residual dependence on M due to the running of the top Yukawa coupling 
as well as the thermal effects included in Abf although it is very mild. 

In Fig. I we plot |?7| as a function of m-e// for M = 10^ GeV. Following Ref. [13,29] 
we consider two different initial conditions for the sneutrino abundance. In one case, one 
assumes that the N population is created by their Yukawa interactions with the thermal 
plasma, and set Yffiz ^ 0) = 0. The other case corresponds to an initial A^ abundance 
equal to the thermal one, Yj;;f{z ^ 0) = Y^^{z ^ 0) §. 

Our results show good agreement with those in Refs. [13-15]. In particular we re- 
produce that for zero initial conditions, r] can take both signs depending on the value of 
meff, thus it is possible to generate the right sign asymmetry with either sign of Im^. For 
thermal initial conditions, on the contrary, r^ > and the right asymmetry can only be 
generated for 1mA > 0. The plot is shown for tan/3 = 30. But as long as tan/3 is not very 
close to one, the dominant dependence on tan/3 arises via Vu as given in Eq. (2.36) and it is 
therefore very mild. For tan/3 ~ 0(1) there is also an additional (very weak) dependence 
due to the associated change in the top Yukawa coupling. 

After conversion by the sphaleron transitions, the final baryon asymmetry is related 



'^The (in)dependence of the final asymmetry on the exact preparation of the initial state has been further 
explored in Ref. [32]. 
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Figure II: B^nrieff regions in which successful soft leptogenesis can be achieved. We take |IniA| = 
10'^ GeV and tan/3 — 30 and different values of M as labeled in the figure. The two panels 
correspond to vanishing initial iV abundance (left) and thermal initial N abundance, (right). 



to the B — L asymmetry by 



Yb 



24 + 4nH 
66 + 13nH 



Yb-l{z -^ oo) 



23 



Yb-l{z 



CXD 



(2.37) 



where nn is the number of Higgs doublets, which is taken to be nn = 2 for the MSSM in 
the second equality. 

This has to be compared with the WMAP measurements that in the ACDM model 
imply [33]: 

Yb = (8.78 ± 0.24) x 10"^^ (2.38) 

We plot in Fig. II the range of parameters B and ?Tie// for which enough asymmetry 
is generated, Yb > 8.54 x 10^^^. We show the ranges for several values of M and for the 
characteristic value of |Imyl| = 1 TeV. 

The figure illustrates our quantification of the known result that independently of the 
N initial distributions, successful soft leptogenesis requires M < 10^ GeV as well as -B ^ ^. 
Next we turn to the effect of flavour on these conclusions. 

Before doing so, let us comment that, as pointed out in Refs. [13-15], these results 
indicate that in soft leptogenesis, the CP asymmetry is maximal when the parameter lie on 
the resonant condition T = 2\B\. In this case, the asymmetry is generated by the decays 
of two nearly mass-degenerate N. It has been recently discussed in Refs. [22] that for 
resonant scenarios, the use of quantum BE [22,34] (QBE) may be relevant. In particular it 
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has been shown that for standard see-saw resonant leptogenesis there are differences with 
the classical treament in the weak washout regime. There, is however, no study on the 
literature of the impact of the use of QBE for the case of soft leptogenesis and to discuss 
those is beyond the scope of this paper. Thus in our work here we study the impact of 
flavour in soft leptogenesis in the context of the classical BE as described above. 

3. Flavour Effects 

In the previous discussion flavour effects have been neglected. This is only justified when 
the process of leptogenesis is completed at temperatures T > 10^^ GeV for which charged 
lepton Yukawa processes are much slower than the processes involving N and than the 
expansion rate of the Universe. 

However, as we have seen, soft leptogenesis is only effective enough for relatively light 
right-handed sneutrino masses M < 10^ GeV. Therefore in the relevant temperature win- 
dow around T ^ M processes mediated by the r and the ^ Yukawa couplings become faster. 
As a consequence, the lepton states produced in the N (and N) decay lose their coherent 
between two subsequent L- violating interactions. So before they can re-scatter in reactions 
involving N and A^ they are projected onto the flavour basis. In this case, the decay rates 
and scattering processes involving the different flavours l^, l^ and anti- flavours l^, l^. have 
to be considered separately and we need to consider the BE for the single lepton-flavour 
asymmetries. 

To account for this effect we need to define the CP flavour asymmetries 



J^7(iV±^afc)-7(iV±^Ofc) 



k _ O-k 



e = ik ^ ^ > (3-1) 



^-f{N±^ak) + 7{N, 



ak 



o-k 



Taking into account that the Yukawa couplings can be chosen to be real the flavoured decay 
rates verify (neglecting the zero temperature masses) 

ak) 

ak) (3.2) 



(3.3) 



with projections K^ 



7(7V±- 


k 


7(7V±- 


-^-ak) = Kiy]^{N^ 
k 




W |2 
^k ,r^. 



so 

e'^(T) ^ e'^ABFiT) = KleABpiT) = Kle{T) (3.4) 

with e(T) given in Eq. (2.25). We notice that because of the assumption of alignment 
between the soft supersymmetry-breaking A terms and the corresponding neutrino Yukawa 
couplings, the only source of CP violation is a flavour independent phase. Therefore, 
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unlike in the case of see-saw leptogenesis induced by N decay [6,7], in this "minimal" 
soft leptogenesis scenario, it is not possible to have non-zero flavour asymmetries with a 
vanishing total CP asymmetry. 

In Ref. [6] the relevant equations including flavour effects associated to the charged- 
lepton Yukawas were derived in the density operator approach. One can define a density 
matrix for the difference of lepton and antileptons such that pkk = ^L,.- As discussed 
in Ref. [6, 7] as long as we are in the regime in which a given set of the charged-lepton 
Yukawa interactions are out of equilibrium, one can restrict the general equation for the 
matrix density p to a subset of equations for the flavour diagonal directions pkk- In the 
transition regimes in which a given Yukawa interaction is approaching equilibrium the off- 
diagonal entries of the density matrix cannot be neglected [7,21]. However, as we will see 
below, for the case of soft leptogenesis, this is never the case. 

There are additional flavour effects associated to the neutrino Yukawa couplings as 
discussed in Ref. [19] such as those arising from processes mediated by N2 and N2 (and 
A^3 and N^). These effects are particularly important in see-saw resonant leptogenesis in 
which right-handed neutrinos of different "generations" are close in mass. In leptogenesis 
with strong hierarchy among the masses of the different generations of right-handed neu- 
trinos/sneutrinos (as we are assuming here) one can neglect the neutrino Yukawa couplings 
in most of the parameter space, because the charged-lepton Yukawa rates are faster at the 
temperatures when the asymmetry is produced. In what follows we will work under this 
assumption and neglect flavour effects associated to the neutrino Yukawas. 

In writing the BE relevant in the regime in which flavours have to be considered, it is 
most appropriate to follow the evolution of lAfe where A^ = -3 — Yl^. — Yl^.^^ = "3 ~ ^L*" • 
This is so because A^ is conserved by sphalerons and by other MSSM interactions. In 
particular, notice that the MSSM processes enforce the equality of fermionic and scalar 
lepton asymmetries of the same flavour. Hence, we can write down the flavoured BE for 



sHz 



dz 



^'iT)i^-^-2y,-Y^M,^,f 






E^fci^TU^ 2 + 



N "N 

(3.5) 



^A, / l^^tot],,fe 

2Yr \ 2 Y'" ' ^22 



j - " \ - ^N 



while the BE for Y^ and Y^ remain the same. 

total 
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In Eq. (3.5) we have defined the flavoured thermal widths 

7^ = i^fe°7^, (3.6) 

l?' = Kh? (3.7) 

The value of A^j^ will depend on which processes are in thermal equilibrium when 
leptogenesis is taking place. For T < (1 + tan^ (3) x lO^GeV where the processes mediated 
by all the three charged lepton (e, /i, r) Yukawa couplings are in equilibrium i.e. they are 
faster than the processes involving N±^ we have [20] 



^ - 40 "M M • (3-8) 



93 


6 


6 


110 


55 


55 


3 


19 


1 


40 


30 


30 


3 


1 


19 


40 


30 


30 



541 


152 


761 


761 


46 


494 


761 


761 



For (1 + tan^ (3) x lO^GeV < T < (1 + tan^ (5) x lO^^GeV where only flavours e + ^ and r 

are distinguishable, we have 

/ _541 152 \ 

(3.9) 

For T > (1 + tan^ (3) x lO^^GeV when all the flavours are indistinguishable i.e. the charged 
lepton Yukawa processes are much slower than the processes involving N±^ we recover the 
unflavoured case where A = —1. 

From the results in Sec. 2 we see that in the relevant temperature window around 
T ^ M and for 1 < tan /3 < 30 we are we are always in the regime when T < (1 + tan^ (3) x 
lO^GeV thus we need to consider flavour effects associated to the three lepton flavours 
separately with A given in Eq. (3.8) 

4. Results 

We parametrize the asymmetry generated by the decay of the singlet sneutrino states in a 
given flavour as 

Ya, (2^00) = -27?, P Y^iT »M) (4.1) 

where e^ is defined in Eq.(3.4). Thus the final total asymmetry can be written as Eq. (2.35) 
where now 

V.fia = Y.^jK^ (4.2) 



In Fig. Ill we plot Irjfia/Tjol as a function of wie// for M = 10^ GeV and for Kf = K\ 



K^ = 3. We label rjQ the corresponding efficiency factor without considering flavour effects. 
As seen in the figure for these values of the flavour projections, Kf, and large m-e// (large 
washout region), flavour effects can make leptogenesis more efficient by up to a factor of 
the order 30. On the contrary flavour effects play no role for small rrieff (small washout). 
This can be easily understood by adding the equations for the three flavour asymmetries. 
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10 10 

nieff (eV) 

Figure III: Efficiency factor |?y/?yo| as a function of nieff for M — 10^ GeV and tan/3 = 30. The 
two curves correspond to vanishing initial N abundance (solid black curve) and thermal initial N 
abundance, (dashed red curve). 



Eq. (3.5). We get an equation which can be written as: 

dYB-L 



sHz- 



dz 



e{T) 



ye? 

, N 



2 ) 7;v - E ^^^- ^° 
/ kj 



2Yr 



eq 



(4.3) 



where we have defined the washout term 



Y, 



Tir , "^tot (5) , r, (6) , n (7) , Y^ (3) , r, (4) , 1 , (2) 



y: 



N 



Y 



N 



l^^totjs) 



„(9) 



Yn (0) 



„(i) 



„(2) 



l^N, 



tot 



+ 2 ^#^^* + 27r^ + 2^7r^ + 27r^ + 27r^ + I 2 + - ^,, 



722 (4.4) 



which can be directly compared with the unflavoured equation Eq. (2.34). We see that if 
we define Pj = Y/^^./Yb-l, Eq. (4.3) is equivalent to Eq. (2.34) with 



W 



-W X Y^^H^jPi 



(4.5) 



Thus flavour effects are unimportant when the W term in Eq. (4.3) is much smaller than 
the source term which happens when me// is small enough (small washout regime). 

We have verified that the equally distributed flavour composition Kf = K2 = K^ = 
1/3 (so all flavour are in the same washout regime) gives an almost maximum flavour effect 
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Figure IV: B,meff regions in wliicli successful soft leptogenesis can be achieved wtien flavour 
effects are included with Kf = K^ = K^ = 1/3 We take |Imyl| = 10^ GeV and tan/3 = 30 and 
different values of M as labeled in the figure. The dashed contours are the corresponding ones when 
flavour effects are not included. The two panels correspond to vanishing initial N abundance (left) 
and thermal initial N abundance, (right). 



for meff < 10^^ eV. Conversely for meff ^ 10^^ eV values, flavour effects lead to larger 
B — L for more asymmetric flavour compositions. In this case, the "optimum" flavour 
projection strongly depends on the value of 'meff- 

We plot in Fig. IV the range of parameters B and rueff for which enough asymmetry 
is generated, Yg > 8.54 x 10~^^ for the the equally distributed flavour composition K^ = 
K2 = K^ = 1/3. We show the ranges for several values of M and for the characteristic 
value of |ImA| = 1 TeV. The dashed contours are the corresponding ones when flavour 
effects are not included. The figure illustrates to what extent flavour effects can affect the 
ranges of B and M for which successful soft leptogenesis can be achieved. This is more 
quantitatively displayed in Fig. V where we plot the asymmetry that can be achieved for a 
give value of B (or M) maximized with respect to w-e// and M (or B) when flavour effects 
are included (for Kf = K2 = K^ = 1/3) compared to the corresponding one when they are 
neglected. From the figure we read that successful soft-leptogenesis with (without) flavour 
effects considered requires B < 8 x 10"^ (3 x 10"^) TeV and M < xlO^ (4 x 10*^) GeV 
for vanishing initial N abundance and i? < 1.5 x 10^^ (3 x 10^'^) TeV and M < 3 x 10^ 
(2 X 10^) GeV for thermal initial N abundance. 

In summary, in this work we have studied the impact of flavour in soft leptogenesis. We 
have quantified to what extent flavour effects, which must be accounted for in the relevant 
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B/1TeV 



M (GeV) 



Figure V: Maximum baryon asymmetry be achieved as a function of B (left) and M (right). The 
sohd (dashed) hncs are for no flavour effects and vanishing (thermal) initial N abundance. The 
dotted (dash-dotted) lines are the corresponding asymmetries after including flavour effects with 
Ki = K2 — K^ — 1/3. The horizontal line correspond to the la WMAP measurements in the 
ACDM model Eq. (2.38). 

sneutrino mass range, can affect the region of parameters in which successful leptogenesis 
induced by CP violation in the right-handed sneutrino mixing is possible. We find that 
for decays which occur in the intermediate to strong washout regimes for all flavours, the 
produced total B — L asymmetry can be up to a factor 0(30) larger than the one predicted 
with flavour effects being neglected. This enhancement, permits slightly larger values of 
the required lepton violating soft bilinear term B. 
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